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Recently, it has been shown that two-dimensional bosonic symmetry-protected topological(SPT) 
phases with on-site unitary symmetry G can be completely classified by the group cohomology class 
H 3 (G, U(l)). Later, group super-cohomology class was proposed as a partial classification for SPT 
phases of interacting fermions. In this work, we revisit this problem based on the mathematical 
framework of G-extension of unitary braided tensor category(UBTC) theory. We first reproduce 
the partial classifications given by group super-cohomology, then we show that with an additional 
H 1 (G, Z 2 ) structure, a complete classification of SPT phases for two-dimensional interacting fermion 
systems for a total symmetry group G x can be achieved. We also discuss the classification of 
interacting fermionic SPT phases protected by time reversal symmetry. 


Introduction — Due to the Pauli exclusion princi¬ 
ple, non-interacting fermions have rich structures in their 
ground state wavefunction and a complete classification 
of symmetry protected topological(SPT) phases of free 
fermions has been achieved by using ideas like Anderson 
localization jlj and K theory Q. It is then of no sur¬ 
prise that the experimental observation of these nontriv¬ 
ial states is relatively easy and a variety of materials have 
been discovered to be time-reversal-invariant topological 
insulators/superconductors On the other hand, re¬ 

cent studies show that SPT phases also exist in interact¬ 
ing boson systems and can be systematically classified by 
(generalized) group cohomology theory [7|-[9j. Neverthe¬ 
less, bosonic SPT phases remain quite elusive since their 
occurrence requires strong interactions (weakly interact¬ 
ing bosons tend to condense and become a superfluid). 

Despite the remarkably successful classifications of 
non-interacting band Hamiltonians, the non-perturbative 
effects of interactions in fermionic SPTs still remain an 
important theoretical question. In many cases, the free- 
fermion classifications are shown to be stable against 
interaction effects, e.g. the Z 2 classification of time- 
reversal-invariant topological insulators [Tp]. A break¬ 


through by Fidkowski and Kitaev HU demonstrated 


that in one-dimensional fermionic systems with time- 
reversal symmetry T 2 = 1, the non-interacting Z clas¬ 
sification breaks down to Zs when strong interactions 
are present 11, 12j]. The result was then generalized to 
two dimensions with an on-site Z 2 symmetry [Tdl-ITsl and 
three-dimensional time-reversal-invariant topological su¬ 
perconductors mu. Recently several general classi¬ 


fication schemes have been proposed for fermionic SPT 
phases, including group super-cohomology [24[, cobor- 
dism 251 and invertible topological field theories J26|. 

In this work we pursue an alternative route to clas¬ 
sify fermionic SPT(fSPT) phases in 2D. Following previ¬ 
ous works on bosonic SPT phases, to characterize fSPT 
phases we introduce extrinsic defects carrying symmetry 


fluxes into the fSPT state. The classification is obtained 
by studying the topological properties of the defects, such 
as their fusion rules and braiding statistics. Similar ideas 
have proven to be quite successful in classifying bosonic 
SPT phases in 2D [ 27 M 30 I ] . The mathematical objects 
that classify 2D fSPT phases can be summarized as three 
group cohomologies of the symmetry group: U 1 (G, Z 2 ), 
BH 2 {G, Z 2 ) and H 3 {G, U(l)). H\G,Z 2 ) which one- 
to-one corresponds to the Z 2 subgroups of G, classifies 
fSPT phases with Majorana edge modes. BH 2 (G, Z 2 ) 
has been derived in the group super-cohomology classifi¬ 
cation 24j and we clarify its physical meaning as the pro¬ 


jective representations carried by fermion-parity fluxes. 
H 3 (G, U(l)) is the well-known classification of bosonic 
SPT phases 0. 

Generalities — First of all, let us clarify the meaning 
of the symmetry group G in a fermionic system. A fun¬ 
damental symmetry of fermionic systems that can never 
be broken is the conservation of total fermion parity, de¬ 
noted by Z^. In addition to this symmetry, we assume 
the system has an on-site symmetry group G. The total 
symmetry group of the system is actually x G [48j . 
For most of our paper we assume G is unitary and finite. 
We will consider anti-unitary time-reversal symmetry in 
the end. 

Our approach to the problem is based on the algebraic 
theory of two-dimensional gapped quantum phases [311. 
Given a gapped phase in two dimensions, one can clas¬ 
sify the low-energy localized quasiparticle excitations into 
superselection sectors (topological charges). Different 
topological charges can not be transformed into each 
other by applying local operators. Fusion and braiding 
of the topological charges are described by the math¬ 
ematical framework of unitary braided tensor category 
(UBTC) [Hj . In a fermionic system, states with even and 
odd numbers of fermions belong to different superselec¬ 
tion sectors. Therefore we model a gapped fermion SPT 
phase abstractly by the (premodular)UBTC C/ = {/, if}, 
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where ip represents a single fermionic excitation. The fu¬ 
sion rule is obviously ip x if = I. 

To classify fSPT phases with a unitary symmetry 
G, we exploit the idea of weakly gauging the symme¬ 
try [17, HI, El]- To be more concrete, we introduce 
into the system non-dynamical point-like defects carrying 
symmetry fluxes. In other words, these defects are en¬ 
tirely extrinsic. We then enlarge our algebraic theory of 
quasiparticle excitations to include these defects. Defects 
by construction carry group labels, but they can also have 
their own topological charge labels. We therefore collect 
all defects labeled by the same group element g £ G into 
a g-sector C g , and define the so-called G-extension Cq as 
Cg = © g eG^g- Notice that the g = 1 sector C\ is just 
the original theory C /. Similar to anyon models, the most 
fundamental property of defects is their fusion rules, i.e. 
how the topological charges are combined. But because 
defects also carry group labels, their fusion rules must be 
G-graded, namely for a g £ C g , b h £ Ch, we have 

a g X &h = 5Z"^ag6 h C g h ’ (1) 

C g h 

where A^ E )[ h are non-negative integers indicating the 
number of ways defects a g and &h can combine to produce 
charge c g h- 

In order to completely define the G-extension, we need 
to study more subtle structures, such as the associativity 
of defect fusion and braiding transformations of defects. 
These are captured in a mathematical formalism called 
G-crossed braided fusion category, and we refer the read¬ 
ers to Ref. [29| for a complete treatment. In our case, 
we follow a more physically intuitive argument to avoid 
solving complicated algebraic equations. In particular, 
we take advantage of a special “invertibility” property of 
SPT phases: for each SPT state, there is a unique “con¬ 
jugate” state such that by stacking them up one obtain 
the trivial state. Furthermore, given two SPT states, one 
can stack them together to get another SPT state, which 
is defined as their sum [49(. In other words, there is a 
natural Abelian group structure on SPT states. We will 
denote the Abelian group of fSPT phases with a given 
symmetry group G by Zl. 

Classifying fusion rules - The G-grading structure 
of the fusion rules of defects has a profound conse¬ 
quence : one can show that all sectors C g have the 
same total quantum dimensions: T > g = T> 2 = 2 where 
V g = £ ageCg d 2 as [H, E3]. In our case, V 2 = 2 leaves 
us with only two possibilities: (a) There are two Abelian 
defects in C g and they differ by fusing with ip. We denote 
them by <r g . (b) There is a single non-Abelian defect in 
C g with quantum dimension \f2. We denote it by o g . 

First we show that the possible non-Abelian fusion 
rules have one-to-one correspondence with H 1 (G, Z 2 ), i.e. 
group homomorphisms from G to Z 2 . Assume for both 
g, h £ G the defects are non-Abelian. To be able to 


construct the fusion outcome of cr g x < 7 h £ C g h, we im¬ 
mediately see that the defects in the gh sector must be 
Abelian just to match the quantum dimension. There are 
still three possibilities: cr g x <Th = 2 er gh ,( 7 g x ay, = 2 <r gh 
and d g x ah = cr g h + cr gh . The former two are impossi¬ 
ble for the following reason: Assuming a g x Oh = 2cr gh . 
Using the symmetry of the fusion coefficients, we must 
have CT gh x er g -i = 2<7h- The left-hand side has di¬ 
mension \/2 while the right-hand already has dimension 
2-\/2, which is clearly impossible. So we conclude that 
t j g xa h = a+ h + a gh . 

On the other hand, if the g sector has a non-Abelian 
defect but the h sector has Abelian ones, the only avail¬ 
able fusion rule is u g x cr^ = a g h, implying that the 
defect in the gh sector is also non-Abelian. Obviously, if 
both g and h sectors have only Abelian defects, the gh 
sector only has Abelian ones, too. 

What we have just established is that whether the g- 
sector is non-Abelian or not gives a homomorphism from 
G to Z 2 . The inverse statement is quite obvious. Given 
any such homomorphism, we can write down fusion rules 
accordingly. Physically, a non-Abelian defect with quan¬ 
tum dimension d, = \/2 is associated with an odd num¬ 
ber of Majorana zero modes localized at the defect which 
implies a topological degeneracy when there are multiple 
defects. 

We also need to determine the fusion rules of the 
Abelian sectors G e consisting of all Z 2 -even group ele¬ 
ments. For any g, h £ G e , we need to specify whether 
u g x (7^ is < 7 g h or u gh = ip x er gh . We can generally write 

<xa+=<^ h) =r (s ' h) x< h , ( 2 ) 

where w(g, h) = ±1 and n{ g, h) = 1 ~ 1 Ug’ h ) . Since fusion 
must be associative, comparing (<r g x < 7 © x and cr g x 

( a h x CT k) y ields 

w(g, h)w(gh, k) = u( g, hk)w(h, k). (3) 

Formally, this means that w(g, h) is a Z 2 -valued 2- 
cocycle: u £ Z 2 {G ei T, 2 ). However, we must remem¬ 
ber that the definition of <Tg is completely arbitrary and 
has no physical meaning. One can always swap a + and 
<Tg for each individual g, or more explicitly, we redefine 
<7g = tp ms x <jg . In terms of <7 g , we find 

“(g' h )- (UyJ “(g' h )- (4) 

w g ,hS related by such redefinitions should be considered 
equivalent. The equivalence classes of 2-cocycles [w] are 
classified by the second group cohomology H 2 (G e , Z 2 ). 

Classifying fSPT phases — As we mentioned before, 
to get a complete classification of fSPT phases we need to 
understand the algebraic data of the defects which satisfy 
a set of consistency conditions j29|. Given a particular 
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fusion rule of defects, there may be more than one dis¬ 
tinct set of algebraic data. On the other hand, for certain 
fusion rules it is possible that the consistency conditions 
do not allow any solutions, in which case the fusion rules 
do not correspond to any two-dimensional fSPT phases, 
i.e. there are obstructions @0. 

Let us start from the cases where all the defects are 
Abelian, i.e. we choose the trivial homomorphism from 
G to Z 2 . Such fSPTs will be referred to as Abelian 
fSPTs. The fusion rules of such extensions correspond 
to 2-cocycles u> £ H 2 (G,Z 2 ). We can show that for the 
obstruction to vanish, the following 4-cocycle must be 
trivial (i.e. it is a 3-coboundary): 


0(g,h,k,l) = (-l) 


n(g,h)n(k,l) 


( 5 ) 


We notice that Eq. ([5]) agrees exactly with the result 
of group super-cohomology 241 50]. Following Ref. 24] 
we denote the obstruction-free subgroup of H 2 {G, Z 2 ) by 
BH 2 {G , Z 2 ), and the group of all Abelian fSPT states as 


&+■ 

We briefly sketch the derivation of 0(See supplemen¬ 
tary material for a detailed derivation). The central 
quantity responsible for the obstruction is the F symbols 
for the associativity of defect fusion, defined diagram- 
matic.ally as 



They must satisfy consistency conditions known as the 
pentagon equations. We can determine the structure 
of E-symbols by bootstrapping from the pentagon equa¬ 
tions between if and the defects, and obtain the following 
general parametrization of defect F symbols: 


=Kg,h,k)A” 


(h,k) 


( 7 ) 


Here A = ± labels the two defects in the same sector, 
and v is a U(l) 3-cochain to be determined. Plugging 
0 into the general pentagon equation we get dv = O, 
which implies that O belongs to the trivial cohomology 
class in H 4 (G , U(l)), thus the obstruction vanishing con¬ 
dition. Once the obstruction vanishes, different solutions 
of v are given by 3-cocycles in H 3 (G, U(l)). In fact, we 
can see from the derivation that U 3 (G,U(1)) is a nor¬ 
mal subgroup of with BH 2 (G , Z 2 ) being the quotient 
group. 

We can actually further constrain the group structure 
of Sf, by noticing that the square of the F symbols 0 is 
v 2 , and dzA = 0, i.e. v 2 £ Z 3 (G, U(l)). This implies that 
“adding up” two Abelian fSPTs given by the same class 
in H 2 (G , Z 2 ) results in a bosonic SPT phase labeled by 

v 2 niiii. 


Next we consider the non-Abelian fSPT phases. First 
we show that given a nontrivial Z 2 homomorphism of G 
there exists at least one non-Abelian fSPT, by explicitly 
constructing the extension. The fusion rules of the Z 2 - 
odd G-sectors are fixed by the homomorphism, and we 
choose a trivial 2-cocycle in H 2 {G e , Z 2 ), i.e. w gi h = 1 for 
the fusion rules of the Z 2 -even Abelian sectors. We define 
a map ip from the topological charges of this extension 
to those of the familiar Ising UBTC, which has three 
topological charges {I, if, a}: 

ip(a+) = I,p(a g ) =ii,ip((T h ) = d. ( 8 ) 

All the algebraic data follow from this map and the data 
of the Ising category. Wc will refer to the corresponding 
fSPT state as the root non-Abelian fSPT state. We now 
argue other non-Abelian fSPT phases with the same Z 2 
homomorphism can all be generated from the root phase. 
We make use of the fact that a non-Abelian defect must 
localize an odd number of Majorana zero modes. Con¬ 
sider two non-Abelian fSPT states fSPTi and fSPT 2 cor¬ 
responding to the same Z 2 homomorphism of G. Denote 
their sum fSPT 3 = fSPTi +fSPT 2 . Suppose we create a 
g-defect. If g is a Z 2 -even element, the defect is already 
Abelian both in fSPTi and fSPT 2 , so is in fSPT 3 . if g is a 
Z 2 -odd element, because fSPTi and fSPT 2 have the same 
Z 2 homomorphisms the defect localizes an even number 
of Majorana zero modes altogether, and can only be an 
Abelian one. Therefore all defects in fSPT 3 are Abelian. 
It immediately follows that any non-Abelian fSPT state 
is equivalent to the sum of the “root” state of the same 
Z 2 homomorphism and an Abelian fSPT state. This pro¬ 
vides a complete classification of the non-Abelian fSPT 
phases. In addition, our argument implies that Sf + is a 
normal subgroup of with H 3 (G, Z 2 ) being the quotient 
group. 

Let us comment on the group structure of non-Abelian 
fSPT phases. From our discussion we see that the root 
non-Abelian fSPT is essentially a Z 2 fSPT. We also no¬ 
tice that given a homomorphism from G to Z 2 , we can 
define a Z 2 2-cocycle on G by pulling back the nontriv¬ 
ial 2-cocycle in H 2 { Z 2 ,Z 2 ). We thus conjecture that the 
addition of two root non-Abelian fSPT with the same 
homomorphism yields an Abelian fSPT given by this 2- 
cocycle. 

Gauging the fermion parity — An interesting probe of 
the symmetry-protected states is to consider gauging the 
fermion parity, which results in a (bosonic) topological 
gauge theory. The gauging procedure is well-defined in 
lattice models of fermions. One can choose the Hamil¬ 
tonian such that the low-energy spectra of the resulting 
gauge theory coincide with the original model. 

For fSPT states, it is easy to see that the gauged the¬ 
ory has four anyons {/, e, m, ip}, where m is the Z 2 gauge 
flux, if is the fermion, e = m x if can be considered as the 
(bosonic) Z 2 charge. The topological order is identical to 
that of a Z 2 toric code lattice model. Importantly, the 
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gauged theory preserves the G symmetry and is there¬ 
fore a G-symmetry enriched Z 2 gauge theory. Therefore, 
a symmetric adiabatic path between two fSPT phases 
maps exactly to a symmetric adiabatic path between the 
corresponding Z 2 topological phases, and consequently 
if the two symmetry-enriched Z 2 topological phases are 
distinct, the original fSPTs must be distinct too [5l|. In 
this approach we can consider anti-unitary symmetries, 
or the fermions carrying projective representations of the 
symmetry group. 


Symmetry enrichment in the toric code model can be 
analyzed using the general theory developed in Ref. 0 
(See also jH] for related discussions). First of all, one 
needs to specify the symmetry action on the topological 
charge labels of anyons. It is easy to see that besides a 
trivial action, there is a Z 2 action that permutes the e and 
to particles. These two possible actions on the label set 
form a Z 2 topological symmetry group. Therefore, the 
symmetry action on the label set is specified by a group 
homomorphism from G to Z 2 . It is clear that this is the 
same Z 2 homomorphism that classifies the non-Abelian 
fusion rules of symmetry defects in fSPT states. 


Once the symmetry action p on charge labels is spec¬ 
ified, we can classify patterns of symmetry fractionaliza- 
tion, i.e. anyons carrying projective representations of 
the symmetry group. This is captured by the group co¬ 
homology Hp(G, Z 2 xZ 2 ), where the subscript p indicates 
that G has a nontrivial action on the coefficients. How¬ 
ever, we need to make sure that the 0 particle always 
transforms trivially under the symmetry. This means 
that we can actually restrict the coefficient of the second 
group cohomology to be Z 2 = {/, 0}, and we drop the 
G action p. Therefore we recover the previous classifi¬ 
cation [37[. We notice that this also provides a physi¬ 
cal characterization of H 2 (G, Z 2 ), through local projec¬ 
tive symmetry actions on the Z 2 fermion parity flux in 
a fermionic system, which can be measured in numerical 
simulations 38, 


Examples — In this section we apply the general the¬ 
ory to G = Z n , Z 2 x Z 2 and Z^. 


• G = Z n . We label the group elements of Z„ by 
a = 0, l,...,n — 1 and the group multiplication is writ¬ 
ten additively, i.e. a + b = [a + b] where [a] is a mod 
n. First we have U 1 (Z„,Z 2 ) = Z( n2 ), be. non-Abelian 
fSPT phases exist for even n. 


We then consider the Abelian fSPTs, classified by 
H 2 { Z n ,Z 2 ) = Z( n 2 ). Let us focus on the case of even 
n. Since U 4 (Z ra ,U(l)) is trivial, there are no obstruc¬ 
tions. A representative cocycle in the nontrivial coho¬ 
mology class is given by uj(a,b) = e^^ + ^~^ a+b ^. In¬ 
deed, w([§], [^]) = — 1 which is a gauge-invariant quan¬ 
tity since we are considering Z 2 coefficients. We can 
also explicitly find the fermionic 3-cocyclc: v(a, b, c) = 
e ^ a ([ fc ]+[ c ]-[ fc + c D^ yy e no tice that i/(a,b,c ) is the square 
root of the generating U(l) 3-cocycle in H 3 ( Zjv,U(1)) = 


Zat. This implies that two fSPT corresponding to the 
nontrivial 2-cocycle can be stacked to form the gener¬ 
ating bosonic SPT. Therefore, the Abelian fSPT phases 
with Z„ symmetry form a Z 2n group. 

We can actually completely determine the group struc¬ 
ture and the classification is given by: 



n is odd 
n = 0 (mod 4) . 
n = 2 (mod 4) 


(9) 


• G = Z 2 x Z 2 . We can easily see that iL 1 (G, Z 2 ) = Z 2 , 
corresponding to non-Abelian fSPTs protected by any of 
the three Z 2 subgroups. The Abelian fSPTs are classi¬ 
fied by 7L 2 (Z 2 x Z 2 ,Z 2 ) = Zj. Using the cohomology 
mapping J5]) we find that all these cocycles are free of 
obstructions. Therefore, all these cocycles represent in¬ 
trinsically distinct fSPT phases. Among the eight, four 
are essentially the fSPT phases corresponding to the two 
Z 2 subgroups. The other four are new fSPT phases where 
a Z 2 fermion parity flux carries two-dimensional projec¬ 
tive representations of the symmetry group. 

• G = Zj = {1,T}. We now discuss time-reversal sym¬ 
metry. Our method of symmetry extension does not ap¬ 
ply because the time-reversal symmetry is anti-unitary. 
However, we can still consider gauging the fermion parity 
and study the symmetry action in the toric code. Since 
after gauging one has a spin model, the time-reversal 
symmetry operator satisfies T 2 = 1. We can then distin¬ 
guish two cases, where the fermion 0 is a Kramers singlet 
(T 2 = 1) or a doublet (T 2 = — 1). 

First let us specify the symmetry action p on the charge 
labels. There are two possibilities: (a) T does not change 
charge labels at all. (b) T exchanges e and to. Interest¬ 
ingly, in the latter case 0 must have T 2 = ~ 1 [3 |4Cj ■ 
Therefore we immediately see that there is a 2D fSPT 
with T 2 = — 1 fermions, in which the local fermion par¬ 
ity of a 7r vortex changes under the time-reversal op¬ 


eration [20j, |41|, |42|, and there are no other symmetry 
fractionalization classes due to H 2 ( Z^,Z 2 x Z 2 ) = 0. 

Let us consider the symmetry fractionalization class 
of the trivial action on the charge labels, which is clas¬ 
sified by 7L 2 (Z^,Z2 x Z 2 ) = Z 2 x Z 2 . Physically, the 
four classes correspond to four possible ways of assigning 
T 2 = ±1 to the four charges. For the two classes with 0 
being a Kramers doublet, one of the e or m charges has 
to be a Kramers singlet, which means that the n vortex 
is trivial in the fSPT. Therefore they do not correspond 
to any nontrivial fSPT phases. We are then left with one 
nontrivial fractionalization class with 0 being a Kramers 
singlet, and both e and m being Kramers doublets. How¬ 
ever, we know that this class has an obstruction of being 
realized in two dimensions [43, 3] and thus does not 
correspond to a fSPT in 2D. Together with the fact that 
H 3 ( Z£\U(1)) = 0, we conclude that there is only one 
nontrivial 2D fSPT phase with Z^ symmetry, the class 
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Dill topological superconductor. 
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Supplementary material 


REVIEW OF GROUP COHOMOLOGY 


In this section, we provide a brief review of group cohomology for finite groups. Given a finite group G, let M be 
an Abelian group equipped with a G action p : G —> M. which is compatible with group multiplication. In particular, 
for any g £ G and a, b £ M, we have 

Pg(ab) = p g (a)p g (b). (10) 

(We leave the group multiplication symbols implicit.) Such an Abelian group M with G action p is called a G-module. 

Let w(gi,..., g„) £ M be a function of n group elements g 7 £ G for j = 1,..., n. Such a function is called a 
n-cochain and the set of all n-cochains is denoted as C n (G,M). They naturally form a group under multiplication, 

(w ■ w')(gi,..., g„) = w(gi,..., g„)o/(gi, • ■ ■, gn), (11) 

and the identity element is the trivial cochain w(gi,..., g n ) = 1. 

We now define the “coboundary” map d : G n (G , M) —► G n+1 (G, M) acting on cochains to be 

n 

dw(gi, . . . , g„+l) = P gl [w(g 2 , . . . , gn+r)] x W (_1)3 (gl, . . • , gj-1, gjgj+1, gj+1, . . • , gn+l) X W (_1) " (gl, . . . , gn). 

3 =1 

(12) 

One can directly verify that ddw = 1 for any w £ C n (G, M), where 1 is the trivial cochain in C n+2 (G, M). This is 
why d is considered a “boundary operator.” 

With the coboundary map, we next define w £ C n (G,M) to be an n-cocycle if it satisfies the condition dw = 1. 
We denote the set of all n-cocycles by 

Z;(G, M) = ker[d : G n (G, M) -)• G n+1 (G, M)\ = { u £ G”(G, M) \ dw = 1}. (13) 

We also define u £ C n (G,M) to be an n-coboundary if it satisfies the condition u = dp for some (n — l)-cochain 
p £ G n_1 (G, M). We denote the set of all ?r-coboundaries by Also we have 

£”(G, M) = im[d : G^^G, M) G n (G, M)\ = { u £ G n (G, M) \ 3p £ G n_1 (G, M) : w = d/a }. (14) 


Clearly, U”(G, M) C Z"(G, M) C C n (G,M). In fact, G n , Z", and are all groups and the co-boundary maps 
are homomorphisms. It is easy to see that B™(G, M ) is a normal subgroup of Z™(G 7 M). Since d is a boundary map, 
we think of the n-coboundaries as being trivial n-cocycles, and it is natural to consider the quotient group 


h;(g,m) 


z;{g,m) 

B?(G,My 


(15) 


which is called the n-th group cohomology. In other words, Hp(G,M) collects the equivalence classes of n-cocycles 
that only differ by n-coboundaries. 


OBSTRUCTION TO ABELIAN G-CROSSED EXTENSIONS 


In this section we derive the obstruction to a consistent Abelian G-extension. We first briefly review the algebraic 
theory of symmetry defects, known as the G-crossed braided extension of a braided tensor category [29j . For simplicity, 
we assume all defects (as well as anyons in the original theory) are abelian. The collection of all defects is called the 
G-extension: 


C G = @ c s- 
g£G 


(16) 


We will use the diagrammatic formulation (for a review, see Refs, 
includes: 


umi). 


The basic data of the G-extension 





• G-graded fusion rules, i.e. a g x bh = 2 Cgh ^C g b h c gh- 

• F symbols for associativity of fusion. 
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The F symbols can be viewed as changes of bases for the states associated with quasiparticles. 


• G action on labels, defined by maps p g , which acts on the topological charge labels in the following way: 
Ph(%) £ Chgh- 1 - We write ph(®g) as h a g . Diagrammatically, this means that when topological charge lines 
cross, the labels should change accordingly: 


a e &h 


6h 



a s 


(18) 


In our derivation we just need the obvious fact that pi(a s ) = a g . 

• G action on fusion spaces, defined by unitary transformations L r k(a g , b^; c g h): 



(19) 


Importantly, we have the normalization condition U\ = 1. 


• Natural isomorphisms r/ Xk (g, h) on topological charges, which define the projective G actions: 



( 20 ) 


Similarly, r] are normalized: ^(1, h) = r] x (g , 1) = 1. 

• G-crossed R symbols, defined by the following diagrammatic relation: 



These data satisfy a set of coherence conditions. For our purpose, we have to solve the pentagon equations and a 
generalization of hexagon equations, called heptagon equations (see Fig. [lj, to find consistent F and R symbols. 

The pentagon equation leads to the following general parametrization of F symbols: 


F° 


^ 3 =i/(g,h,k)0 2 1 (h, k). 


( 22 ) 








Notice that because we are considering Abelian fusion rules, all labels in the definition ED of F symbol are uniquely 
determined by the three outgoing lines and therefore suppressed here, ^(g, h, k) is a U(l) 3-cochain, and 6 is a Z 2 
2-cocycle. And they should satisfy 


(di/)(g,h,k, 1) = 6>"S’ h (k, 1). (23) 

This implies that the right-hand side must be a 4-coboundary. 



d d d d d d 


FIG. 1: Heptagon equations which ensure consistency between G-crossed braiding and fusion of defects. 


We now use G-crossed heptagon equations to find 6. 

Let us first consider the heptagon equation with the three outgoing lines being < 7 ^, 

R^i 1 p a i^ R^ = R a s 


ip and ip (from left to right). 


(24) 


Since RW = -1, we have 

R^ =-R 7 ^. (25) 

Let us consider the heptagon equation with the three outgoing lines being Cg 1 ,^ 2 and ip (from left to right). 
Because ip € Ci, the action on the vertex is trivial. We have 


= CT h 2 R 7 ^ XCT h 2 ^ p^g 1 A) 2 i/i 


(26) 


Following the general parametrization and the normalization condition of v, we have F ,0 'g 1 ^" T h 2 
F^ 1<T h 2 =6»(g,h). Therefore 


= 1 and 


0(g,h) 


R^^R 7 h 2 ^ 
R 7 g 1 xo-h 2 ,b 


Now we set Ai = A 2 = + and notice <7+ x = ip nl ' s ' h) x a^ h , we arrive at 


0(g,h) 


^—l) n S,h 


R a t^R a U 
R<^ 


(27) 


(28) 


We can further show that R a £^ = ±1 by considering the heptagon equation for inverse braiding with the three 
outgoing lines being ip, ip, a*. This is basically what we need, i.e. 9 is Z 2 -cohomologically equivalent to (—l) n . 
Therefore, the necessary condition for the extension to exist is that the 4-cocycle: 


0(g, h,k, 1) = (— l) n B,h^k,i 


(29) 


is in the trivial cohomology class, i.e. [ 0\ = 0. Once the obstruction vanishes, different solutions of v are related to 
each other by a 3-cocycle. 
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FIG. 2: The Pentagon equation enforces the consistency between different sequences of F moves starting and ending with the 
same fusion trees. 


PROJECTIVE PENTAGON EQUATION 


In this section we give a different derivation of the obstruction mapping, which is closer in spirit to the group 
super-cohomology theory. 

Instead of considering all the defects , we take a representative, e.g. (7+ from each g-sector. To account for the 
fusion rules, we allow the fusion space to be fermionic, i.e. Z^-graded. In our case, the fermion parity of the fusion 
space of cr g and Oh is completely determined by g and h: it is just (— l) n (s. h ). More general constructions of fermionic 
TQFTs have been studied in Refs. [46, 47). 

With this modification, now the F move becomes an operator possibly connecting states with different fermion 
parities. 



= ^(g,h, 


k) gh,k) j.n{ g.h) f-n( h,k) j-n( g,hk) 



(30) 


Here when we write / we really mean the fermionic creation operator while / _1 means the annihilation one. 

Let us substitute the fermionic F moves into the pentagon equation. We will not explicitly write the v factors 
since they just give the standard 4-coboundary, and focus on the fermionic part. First evaluate the upper path of the 
Pentagon equation: 

J7l(glg2g3,g4) yn(glg2,g3) y — ra(g3,g4) y -n(glg2 ,g3g4 ) y™(glg2 ,g3g4 ) y™(gl,g2) y — 7l(g2 ,g3g4 ) y-n(gl ,g2g3g4) 

7 A J 3 7 2 J A 7 4 71 7 5 7 4 . 

_ y«(glg2g3,g4) rn( glg2,g3) y — 7l( g3,g4) y«( gl.g2) y-n(g2 ,g3g4 ) y -n(gl ,g2g3g4) ' 

— JA J 3 J 2 /1 J 5 JA 


Then the lower path gives: 


f n(glg2,g3) yn(gl,g 2 ) y-n(g2,g3) y-n(gl,g2g3) v y»(glg2g3,g4) yn(gl,g2g3) y-n(g2g3,g4) f~n( gl,g2g3g4) 

73 71 76 73 A J 4 J 3 J 5 J 4 

v y7l(g2g3,g4) yn(g2,g3) y-™(g3,g4) y-n(g2,g3g4) 

a 75 76 72 75 


(32) 


We move the first j*( glg2g3 - g4) t 0 the left most and then / 4 s ( gl,g2g3g4 ' to the right most and the resulting fermionic 
sign is the following: 


(- 1 )' 


"(glg2g3,g4)[7l(glg2.g3) + n(gl.g2)-n(g2.g3)-n(gl.g2g3)]^_t j- n (gl.g2g3g4)[7l(g2g3,g4)+n(g2,g3)-n(g3.g4)-7l(g2.g3g4)] 


(33) 


From the 2-cocyle condition 


w(gl, g2g3)w(g2, g 3 ) = w(gig 2 , g 3 )w(gl, g 2 ) 


( 34 ) 
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together with oj = (—1)”, we have 

n( gi, g 2 g 3 ) + n( g 2 , g 3 ) = n(gig 2 , g 3 ) + n{ gi, g 2 ) (mod 2). (35) 

Therefore the term (l33l) is identically 1. 

After this manipulation the expression (l32l) can be greatly simplified: 


gig 2 g3,g4) fn(gig 2 ,g3) j?n( gi,g 2 ) ir-n(g 3 ,g 4 ) f-n( g 2 ,g3g4) *-n(gi ,g 2 g3g4) 

./ 4 J 3 J 1 J 2 J 5 J 4 


(36) 


So the upper and the lower paths (l36l) just differs by exchanging y" lgl ' g2 l and "' g3 ' g4 > which results in a sign: 

0(gi,g2,g 3 ,g 4 ) = (- 1 )™^^^). ( 37 ) 

In order to satisfy the pentagon equation, this sign has to be compensated by dv: 

(di/)(g, h, k, 1) = (— i) ra (S’ h )"(k i)_ ( 38 ) 


This is the obstruction-free condition. 


SYMMETRY FRACTIONALIZATION IN Z 2 TORIC CODE 

We derive the H 2 (G, Z 2 ) classification through the symmetry fractionalization in Z 2 toric code. Our argument here 
is somewhat heuristic and we refer the readers to Ref. [29 ] for more rigorous discussions. 

To be more specific, let us consider the symmetry action on a general quasiparticle state on a sphere/disk without 
loss of generality, and assume that the permutation p is trivial for simplicity. The global symmetry operator R g can 
be decomposed into operators localized on each quasiparticle U g (a). They only form projective representations of G, 
i.e. U g (a)Uh(a) = rj a (g,h)U g i 1 (a), however rj a (g, h) must be consistent with fusion rules: 

Va( g, h)%(g, h) = ?7c(g, h), if N£ b > 0. (39) 

This stems from the fact that R g Rh = i? g h must hold on vaccume state. In particular, we have rj 2 = r]'^ = 1, 77 ^, = 

VeTlm- 

On the other hand, from the associativity of operator products U g UhUk , we have 

7? Q (h,k)77 Q (g,hk) = 77 Q (g,h)77 a (gh,k). (40) 

So naively, one may conclude that gauge-inequivalent classes of 77 are classified by H 2 (G, Z 2 x Z 2 ). However, since 
we are considering symmetry-enriched toric code from gauging a fSPT phase, the if} quasiparticle has to form a linear 
representation, implying rj^ can be chosen to 1 and therefore r/ e = rj m . So the actually classification is just H 2 (G , Z 2 ). 

We also notice that any Abelian phase on anyon a that satisfies a relation like (l39l) must be the braiding phase of 
a (fixed) Abelian anyon with a, which actually holds for any topological phase following from modularity. The fact 
that we have 7? e = rj m means r] e ( r]m ) is the braiding phase of either I or ip with e(m). Physically, we can think of 
U g (a) as taking the cr+ defect around a. If rj e = rjm = —1, the only possiblity is that when er+ fuses with er^, we 
obtain a^ h and also a fermion ip which when taking a full braid around e(m) generates the phases r/ e (r/ m ). Therefore, 
if we write ?7e(g, h) = r} m (g, h) = (— l) n (s> h ) with n( g, h) = 0,1 being the linearized Z 2 2-cocycle, we have 


a+xa+=r (s ’ h V+ h - 


(41) 


